Spin-boson coupling in continuous-time quantum Monte Carlo 
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A vector bosonic field coupled to the electronic spin is treated by means of the continuous-time 
quantum Monte Carlo method. In the Bose Kondo model with a sub-Ohmic density of states 
Pb(w) ~ lu~ 3 , two contributions to the spin susceptibility, the Curie term T _1 and the term T~ s 
due to bosonic fluctuations, are observed separately. By including the fermionic bath, a quantum 
phase transition is identified between the Kondo screened state and the bosonic fluctuating state, at 
which the effective moment and the local Fermi-liquid energy scale vanish. It is demonstrated that 
the energy scale of the bosonic fluctuation is not affected by the existence of the quantum phase 
transition. 

PACS numbers: 75.20.Hr, 71.10.-w 



I. INTRODUCTION 

The continuous-time quantum Monte Carlo (CT- 
QMC) method has been developing since 2005, as a 
numerical tool for correlated fermion systems. [l|, 0] Es- 
pecially, the algorithm based on the expansion around 
the atomic limit (CT-HYB) is highly effective as the 
impurity solver for the dynamical mean-field theory 
(DMFT).041 The method can also be applied to vari- 
ants of Kondo models, where a localized spin interacts 
with the fermionic bath via the exchange coupling. [(| 0] 

There is another class of impurity models which con- 
sist of an additional coupling between a bosonic bath 
and local degrees of freedom. The simplest one is the 
coupling between the electronic charge n/ and a boson 
<f> of the form rif(f). In CT-QMC, arbitrary energy dis- 
persion of the bosonic bath is treatable, and a dynamical 
screening effect has been investigated. [8j This algorithm 
can also be applied to the coupling Sjcj) with 5| being 
^-component of the local spin.Q We may consider more 
complicated interaction including a spin flip scattering of 
the bosonic field, leading to the coupling Sf ■ 4>, where a 
vector bosonic field couples to the electronic spin Sf. 

The coupling Sf<f> appears when the Heiscnbcrg inter- 
action is treated in a "mean-field" theory. The boson </> 
describes an auxiliary field which mediates the effective 
local spin-spin interaction resulting from the intersite in- 
teraction. A "mean-field" treatment of the quantum spin 
glass introduces this interaction with a self-consistent 
equation. 

In the (non-rando m) Hcisenberg model, 
this impurity model describes spin fluctuations around 
the molecular field in magnetic ordered states as well as 
in the paramagnetic state. 14| With a fermionic bath in 
terms of DMFT,[l5[ doping of the spin glass [l6| as well as 
the Heisenberg-Hubbard model can be addressed beyond 
the molecular- field approximation. fl7l - [l9l ] This extended 
DMFT equation has recently been reformulated in the 
dual boson approach. [2^ A similar model has also been 
investigated in the context of a non-magnetic impurity 
embedded in an antiferromagnet . [2l| 



The impurity Hamiltonian consisting of the fermionic 
bath Ofco- and the vector bosonic bath b q £ (£ = x, y, z) 
reads 



H = E e !° n f° + U n tf n fi + E eka k<J 



O-ka 



E ^ b \M + v E(/- a - + + E 

(i) 



where a, 
N 



• = of -|- Ug, and = 
1 / 2 ^2 q b q £ with N being the number of sites. rif a = 

fife, and Sf = J2<j<j> ft^fv' with ^ being the Pauli 
matrix. We introduce the XXZ-type anisotropy for the 
bosonic bath and the spin-boson coupling, u qx = uj qv = 
uj ql _ and g x = g y = g±. 

The bosonic part in H is reminiscent of the spin-boson 
model, which has been investigated in the context of dis- 
sipative systems. [22l |23| Its SU(2) symmetric version is 
referred to as Bose Kondo model. [21( On the other hand, 
the model including a spin exchange coupling with the 
fermionic bath is known as Bose- Fermi Kondo model. [TtJ 
The Hamiltonian |1| describes charge fluctuations as 
well, and may be referred to as Bose-Fermi Anderson 
model. 

The essence of this model is that the fermionic bath 
screens the localized spin, while the bosonic bath stabi- 
lizes the local moment to decouple the fermionic bath. 
This competition, in a certain situation, leads to a quan- 
tum phase transition between the Kondo-screened state 
for small g and a local-moment state with a residual mo- 
ment for large g.[24[ When the SU(2) symmetry for the 
spin is preserved, the local-moment state is governed by 
an intermediate-coupling fixed point. [2l| A critical na- 
ture has been clarified by means of renormalization group 
theory with inclusion of full anisotropy. pET - E^ 

In numerical approach, on the other hand, the case 
of Ising-type coupling has been investigated by numer- 
ical renormalization group (NRG).[lj| In this case, the 
local-moment state is governed by a strong-coupling fixed 
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point. In this paper, we address the isotropic coupling 
with SU(2) symmetry and explore the physical quantities 
arising from the intermediate-coupling fixed point as well 
as the quantum phase transition. 

In the next section, we present an algorithm based on 
CT-QMC to solve the model ([1]). We restrict ourselves 
to U = oo in the simulation, which is related to the t-J 
model and the Hciscnbcrg model in terms of the extended 
DMFT. Numerical results are presented first for a pure 
bosonic system V = (Bose Kondo model) in Sec. Mil By 
including the fermionic bath, the quantum phase transi- 
tion is explored in Sec. IIVI A summary is finally given in 
Sec.lV] 
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FIG. 1: An example of the Monte Carlo configuration of order 
k = 2 and 1 = 2. The dark and light shaded area indicate 
the spin-up and -down states, respectively. The curved lines 
express the bosonic Green function. 



II. SPIN-BOSON COUPLING IN CT-QMC 

We solve the effective impurity model ([J} using the 
hybridization-expansion solver of the CT-QMC. 0, Q In 
this section, we present how to treat the additional 
bosonic field in CT-QMC. 

The bosonic bath coupled to the electronic charge has 
been formulated by Werner and MillisQ. In this method, 
the electron-phonon coupling is eliminated by the so- 
called Lang-Firsov transformation, and it makes the com- 
putation efficient. This manipulation can also be applied 
to the coupling between Sj and bosons. In the pres- 
ence of the spin-flip term, however, the above treatment 
cannot be applied, since the spin operators, SJ, Sj, and 
Sj, do not commute with each other. Only one compo- 
nent can be eliminated among three. We evaluate the 
other two couplings by a stochastic method. Namely, 
hybridization as well as the spin-boson coupling are eval- 
uated via a Monte Carlo sampling. 

Before going to the formulation, we define the prop- 
agators for the fermionic bath (hybridization function) 
and the bosonic bath (effective interaction) as follows: 



A(iw n ) = V 2 G (iw„) 



V 2 



k 



Sy(iv n ) = -g 2 D l0 {iv n ) = -1 ^ 



2u, 
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(3) 



where ui n = (2n + l)7rT and v n = 2nirT are the fermionic 
and bosonic Matsubara frequencies and 7 = z, JL. The 
latter quantity describes the effective interaction —Sf(r)- 
J(t — t')Sj(t') mediated by the bosonic field. 



A. Canonical transformation 

We first eliminate the coupling between Sj and bosons. 
Following Ref. Q , we perform a canonical transformation 

H = e s He- s with S = N' 1 / 2 £ (<W«V)(H* " h ^) S P 
which shifts the z-coordinate of the oscillation to elimi- 
nate the term g z Sj<j) z . The transformed Hamiltonian H 



is given by 

H = l f° n f° + & n n n fl + /Z £kal a a ko 



kcr 



+ E "A h <* + v E(/- a - + 

+ E±(S+^ + Sj<j> + ), 



(4) 



where ± = ((f> x ± i(f> v )/\2. The local parameters are 
renormalizcd to if a = e/ CT — N^ 1 g 2 / (4ui qz ) and 

U = U + N- 1 J2 q gl/^Uqz)- The operators for the local 
electron are transformed to 

U = e~° A / 2 U, fl = *° A/% ft> Sf = e ±A Sf, (5) 

where Sf = Sf ± iS} and A = N' 1 ' 2 E 5 fe/^)(6^ - 

b qz ). In Eq. ([5]), the factor e A ^ 2 is associated with the 
change in the quantum number of Sj. 



B. Partition function 

With the transformed Hamiltonian H, we expand the 
partition function Z with respect to V and <7j_ as follows: 



Z 



1, n J n J J 



(6) 



k=0 1=0 



where the subscript denotes a quantity for V = g± = 0. 
The integrand W{r, fi) describes the contribution of or- 
der V 2k g±. t = (ri, • • • ,T2fc) denotes a set of imaginary 
times at which the hybridization events occur. The inte- 
gral over r is taken for the range /3 > T2k > • • • > Tl > 0. 
(i = (fi±,--- , fJ-21) denotes a set of imaginary times at 
which a spin exchange takes place. The integral over fj, 
is taken in the same manner as t. Figure [TJ shows an ex- 
ample of the configuration. The doubly occupied state is 
excluded in this figure, since we consider the limit U = 00 
in the next subsection. We note that the formulae in this 
subsection are valid also for U < 00. In the simulation, 
the summations over k and I as well as the integrals over 
r and \i are to be evaluated via the importance sampling. 



The weight W(t, /x) is decoupled into four contribu- 
tions: 



W{T,H) = W loc (T,fJ,)W hyh (T)W ± (fJ.)W z (T,fJ,). 



They denote, from left to right, thermal averages over 
the /-state, the fermionic bath, the xy-component of the 
bosonic bath, and the z-component. The first two are 
the contributions from the Anderson model: Q W\ oc is the 
trace over the /-states with the renormalized parameters, 
if a- and U , and Whyb is the trace over the fermionic bath, 
which is expressed by the determinant of k x k matrix 
consisting of A(r) in Eq. (0. 

W±(fJi) describes the contribution from the g^-terms 
and is given by 



gf 



'M-"^(Ml)>0, 



(8) 
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FIG. 2: Update processes necessary to evaluate the spin- 
boson coupling. 



where (ffl denotes either <j) + or </>~ . The numbers of 
(j) + and </>~ must be the same since we have the re- 
lation (^> ± (/Xj)^) ± (/Xj))o = 0. The thermal average in 
W± is decomposed by Wick's theorem, which is repre- 
sented by the permanent of I x I matrix consisting of 

D±o(fi>i — (J-j) = -(T T (/> + (fii)(t>~(Hj))o-Wk However, since 
there is no efficient algorithm to compute the permanent, 
we evaluate it by stochastic sampling.^ Namely, we ex- 
press W± as 



(9) 



with a denoting one of terms in the permanent, and the 
summation is to be evaluated stochastically. In Fig. [TJ 
the configuration a is represented by curved lines. 

W z (t,h) consists of the phase factors in Eq. (|SJ) and 
is represented as 



W z (r,ti) 



D s 2m A(t 2m ) 



■ e SlA(tl) )o, 



(10) 



where {U} is composed of r and \x in ascending order 
and m = k + I. A{t) = N' 1 / 2 J2 q (g z /uj qz )(e^% z - 
e~ u " zt b qz ). The factor s % takes a/2 for /t, -a/2 for 
/j , and ±1 for . Using the condition J2i s i = 0; the 
thermal average can be evaluated analytically to giveQ 



W z (r,fi) = cxp 



J2 SiSjKfe-U) , (11) 

2m>j>i>l 

K ^ = -jjH 4-[ B K». ^ - S K*> °)]- ( 12 ) 

where r) = cosh[(/3/2 - r)a;]/sinh(/3w/2). 

The dynamics of the bosonic bath is described by 
the function J 7 (iv n ) defined in Eq. ([3]). It is there- 
fore convenient to express the summations over q in the 
above formulae in terms of JJ 7 (ii/ n ). The renormalized 
parameters are rewritten as e/ CT = e/ CT — JT z (0)/8 and 



U = U + J x (0)/4. On the other hand, if(r) in Eq. (JH 
is expressed as 



K(t) = J,(0) 



rW - t) 

2/3 



-L^'W — ^ — ■ ( 13 ) 



C. Monte Carlo procedure 

We perform stochastic samplings of r and /x in Eq. ([6]) 
and a in Eq. @ . They respectively correspond to the V- 
cxpansion, g^-expansion and the Wick's theorem for the 
bosonic field. Since the Hamiltonian with V — g±_ = 
conserves the quantum number of Sp we can treat W\ oc 
by the "segment picture" of CT-HYB.0, Q Hence for 
the U-expansion, the update procedure in the Anderson 
model can be used.Q Hereafter, we consider the limit 
U = oo, which can be implemented by excluding the 
doubly occupied state in the configuration. 

In addition to the updates in CT-HYB, we perform the 
following updates: 



(a) Insertion/Removal of Sf(/j 



t)Sf (/i) on t-state. 



(b) Insertion/Removal of Sj (// + l)S~f(n) on U-state. 

(c) Change of the configuration a. 

(d) Replacing S+(^) and Sj(m) with f\{^+l')fM) 
and fl(fJ>i + d)ff(fii), and vice versa. 

These updates are expressed diagrammatically in Fig. [2] 
The updates (a) and (b) change the expansion order of 
g± by 2. In (c), we exchange two links of the bosonic 
Green functions. The ergodicity is in principle satisfied 
only by (a)-(c). However, a part of the configuration 
may freeze in practice, when the expansion orders for g± 
and V arc considerably different from each other, say, 
when g± is much smaller than V. The freezing happens 
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because the spin operators are separated with each other 
by hybridization operators being located in between. In 
such a situation, the spin operators cannot be removed by 
the updates (a) and (b) . This problem can be resolved by 
introducing the update (d), which replaces two spin op- 
erators separated in the time ordering by single-particle 
operators. [For example, this update is important in the 
parameter range < g < 0.12 in Fig. [5J] 

We first consider the update (a). In the insertion pro- 
cess, we choose two imaginary times randomly in the 
same way as the "segment algorithm" [HQ: /Lt is first 
chosen from the range [0 : /3) and then the length I is 
chosen in the range (0 : ^ max ) so that the operator Sj 
does not pass the next operators. In the removal process, 
we choose one pair from (/ + 1) pairs of the spin opera- 
tors which are connected by the bosonic line, and try the 
update if it is allowed, i.e., if no operator exists between 
them. From the detailed balance condition, the update 
probability R is given by 



R{H ^ 



J±{-l)W loc (T,n+)W z (T,fi+) 



l + l 2 W 1oc (t,») W z (t,h) ' 

(14) 

where fi and fi + denote the configurations of order gf_ 

and g± l+1 \ respectively. The expression for the update 
(b) is given in a similar manner. 

The update probability for (c) comes only from W±. 
Suppose that (/Zj,/4) and (fMj,^) denote the pairs of 
imaginary times connected by the bosonic Green function 
in the original configuration a, the update probability R 
of exchanging the links is given by 



R(a 



(15) 



Finally, we consider the update (d). We first choose a 
pair of spin operators connected by the bosonic line, as 
in the removal process in the update (a). They are to 
be replaced by fi and /f, respectively. Simultaneously, 
the operator ft (/[) is placed next to /j_ (ff). Here, the 
length of the empty state I (£') is chosen from the range 
up to £ max (£' max ) so that the resultant configuration is 
allowed. In the opposite process, we choose the operators 
fl and ff from (k± + 1) and (fc-j- + 1) randomly, where 
k G denotes the hybridization-expansion order for spin tr. 
The update probability R is given by 



R(r,p 



T++, fi- 



ll £' 



W hyh (T + +) W loc (r++ , p-) W z (r++,fi-) 



W hyh (r) W loc {r,») W z {t,h) 



, (16) 



where t ++ and fi denote the new configuration of order 
V 2(k+2) g 2Q-i)_ 

We have confirmed, in the simulation, that all the up- 
date probabilities presented above are always positive 
and therefore, the simulation does not suffer from the 
sign problem. 



D. Spin susceptibility 



We define the spin susceptibilities by Xzz(t) = 
(S z f (T)S z f ) and x+-(r) = (S^(t)SJ)/2. In the isotropic 

system, we have Xzzij) = X-\ ( T )- We can evaluate 

Xzz (t) from the configuration of the /-operators as in the 

"segment algorithm" [2, ■ On the other hand, \-\ ( r ) 

can be evaluated by 



X+-(r) 




(17) 



MC 



where fa and denote the imaginary times for Sj and 
which are connected by the bosonic Green function, 
and MC means average over Monte Carlo configuration. 
The function 6(r,fi) is defined by 



6(t-h) (a*>0) 
&(r-ii-P) (ji<0) 



(18) 



and xa ( T ) i s sampled in the range < r < /3. The 

end points are evaluated accurately from the occupa- 
tion number using the relations x-\ (+0) = ( n /t)/2 an d 

X+-ifi - 0) = {n fi )/2. Eq. ([T7J) follows from the fact 
that J± describes the retarded interaction between the 
local spin so that it may be regarded as the source field 
for the susceptibility. 

The susceptibilities Xzz(t) and x-\ ( T ) can a ls° be 

computed using the matrix M G which is kept in the 
simulation to evaluate the determinant in Wh y b ■ Al- 
though this way is not efficient compared to the method 
presented above, we can use it for a check of the algo- 
rithm and a code. Another consistency check is Xzz(j) = 

X-\ (t) in isotropic parameters, since this condition is 

not trivial in the present algorithm, which treats g z and 
g± in different ways. We have confirmed that our results 
satisfy this condition. 



III. PURE BOSONIC SYSTEM 

In this section, we present numerical results for the 
pure bosonic system, V — 0. In our simulation, the 
charge fluctuation is absent because of U = oo so that 
the local electron is reduced to a localized spin S. Hence, 
the corresponding Hamiltonian is written as 



H 



BK 



gS ■ <(>. 



(19) 



This model is referred to as Bose Kondo model, which has 
been investigated in the context of an impurity embedded 
in an antiferromagnet.[2l],[25[ We restrict ourselves to the 
isotropic system, uj qz = ui q ± = uj q and g z = g± = g. 

The bosonic bath is characterized by the density of 
states pb(w) = N~ lJ £2 q S(uj — u> q ). We use a function 
Pb(w) oc uj s with a cut-off energy uj c . The sum- rule of 



T 



g 



FIG. 3: (Color online) Temperature dependences of the static 
susceptibility x(0) in the pure bosonic system with s = 0.2 
(solid lines). The regular part Xrc g (0) defined in Eq. (1211) is 
also plotted (dashed lines). 

the density of states, (kop^ui) = 1, determines the 
factor to yield the explicit form 

p B (w) = (s + l)w s w- s -^(w c - w). (20) 

We take lo c = 1 as the unit of energy. 

According to the renormalization group analysis. (2li 
HE, H3, HE] this model has an intermediate fixed point 
g = g* for < s < 1. At this fixed point, the suscepti- 
bility shows the long-time behavior x( r ) ~ T > which 
indicates the static susceptibility of the form x ~ T~ a . 
For s > 1, on the other hand, the bosonic coupling is 
irrelevant. We investigate properties due to the interme- 
diate fixed point restricting ourselves to s = 0.2. 

Fig. [3] shows temperature dependences of the static 
spin susceptibility for s = 0.2. We see the Curie behav- 
ior x T at low temperature indicating existence of 
a residual moment. To see the singularity T~ a originat- 
ing from the intermediate fixed point, we define a regular 
part of the susceptibility, Xreg(z), by an analytical con- 
tinuation of x( z = i v n) with v n > 0. Using Xreg(z), the 
susceptibility is written as 

X{w n ) = S n0 M/4T + XregOn)- (21) 

We note that the effective moment M may depend on 
temperature. The full-moment corresponds to M = 1. 
We evaluate Xreg(0) by an extrapolation from x(^i), 
x(«^2) and xiivs) with a quadratic function. We have 
confirmed that the choice of the functional form in 
the extrapolation does not affect the low-temperature 
behavior. [32j The result is shown in Fig. |3l We clearly 
see the power-law behavior Xreg(0) oc T~ s at low temper- 
atures. Consequently, the low-temperature static suscep- 
tibility can be expressed in terms of two diverging terms 

X(0) ~ Mo^T+I/^'CTb) 1 -], (22) 



FIG. 4: (Color online) The energy scale Tb of the bosonic 
fluctuation in the pure bosonic system with s = 0.2 (denoted 
by BK). The result for the Bose-Fermi Anderson model is also 
plotted (denoted by BFA, see Fig. [5] for parameters). The 
solid line is the function Tb oc g°" fitted to the BK data. 

where Mq = liniT-^o M is the residual moment and we 
have introduced a characteristic energy scale Tb of the 
bosonic fluctuation. It turns out that Tb exhibits a 
power-law behavior as shown in Fig. 01 The exponent 
a defined by Tb oc g a is obtained as a ~ 2.52 with the 
error 0.01. The residual moment Mo weakly depends on 
g. [The figure is presented in the next section (Fig. [6j to- 
gether with results for the Bose-Fermi Anderson model.] 

IV. FERMIONIC AND BOSONIC BATHS 

We proceed to system with the bosonic and fcrmionic 
baths. For the bosonic bath, we fix s = 0.2 as in the 
previous section. For the fcrmionic density of states, 
Pf(cj) = N^ 1 J2k — e fc)> on * ne other hand, we use a 
rectangular model with a cut-off energy D 

p F (u) = (l/2D)0(D-\u\). (23) 

We vary g z = g± = g fixing V 2 — 0.1, e/ CT = —0.2, 
U = oo and u) c = D = 1. The Kondo temperature Tk is 
estimated to be Tk ~ 0.1 for g = 0. 

According to the renormalization grou p analysis for 
the Bose-Fermi Kondo model, [2^, l27l |28| there arc two 
fixed points which are stable within SU(2) symmetry in 
(J,g) plane, where J = V 2 /\ef a \ in the present model. 
One is the fermionic fixed point with (J,g) = (oo, 0), at 
which the bosonic bath is decoupled to form the Kondo 
state. The other is bosonic fluctuating fixed point with 
(J, g) = (0,g*), which has been examined in the previous 
section. A quantum phase transition takes place between 
these two states. At the critical point, the susceptibility 
shows the power-law behavior x ~ T s . 

We show temperature dependence of the spin suscep- 
tibility in Fig. [S] Difference with the pure bosonic sys- 
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T 



FIG. 5: (Color online) Temperature dependences of the static 
susceptibility x(0) in the Bose-Fermi Anderson model with 
V 2 = 0.1, e ftT = -0.2, U = oo and s = 0.2 (solid lines). The 
dashed lines show the regular part Xreg(0) defined in Eq. ([21) . 

tern is the paramagnetic behavior in the small-g region, 
g < 0.10. This indicates the spin fluctuation in the 
Kondo singlet state given bv(33j 

X (0) = 1/4T F , (24) 

with Tp being the energy scale of low-energy excita- 
tions. The susceptibility increases against g indicating 
a reduction of Tp- To quantify the local Fermi- liquid 
state, we evaluate the renormalization factor z defined 
by z = [1 — ImS/(ia;o)/wo] _1 - The result is plotted in 
Fig. |5] for several values of T. A reduction of z with in- 
creasing g is demonstrated and at g = 0.12 (0.14), it is 
estimated as z < 0.023 (0.008). However, since z is not 
converged around g = 0.12 in this temperature range, 
we cannot identify the quantum critical point from this 
data. 

In the large-g region, g > 0.14, on the other hand, 
X shows the Curie behavior x x T . As in the pure 
bosonic system, we evaluate the regular part Xreg(0) de- 
fined in Eq. (|2ip . which expresses the contribution ex- 
cepting the Curie term. It turns out from Fig. [5] that 
Xrcg(0) shows the power-law behavior T~ s as in Fig. [3] 
A remarkable point is that the energy scale Tb of the 
bosonic fluctuation is not affected by the fermionic bath 
as shown in Fig. Hence, the difference to the pure 
bosonic system in the local-moment regime comes from 
the Curie term. To see this, we evaluate the effective 
moment M by subtracting Xrog from x in Eq. ([2~Tj) . and 
plot it as a function of g in Eig. [SJ It turns out that 
M is strongly suppressed compared to that in the pure 
bosonic system below g ~ 0.20. 

Finally, in order to identify the quantum critical point 
g c , we plot T s x as a function of g for different tempera- 
tures in Fig. [7J The low-temperature expression of x m 
Eqs. (f2"2")l and (j2~4")l indicates that T s \ is independent of 
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FIG. 6: (Color online) The renormalization factor z and the 
effective moment M for the parameters same as in Fig. [5] Re- 
sults at four different temperatures are plotted. The dashed 
line (denoted by BK) is the result for the pure bosonic system. 
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FIG. 7: (Color online) T B \ as a function of g for different 
values of T. From the intersection points, which are indicated 
by arrows, the critical point is determined. The inset shows 
extrapolation of the intersection point to T = 0. 



temperature at g = g c . while T s \ = for g < g c and 
T s \ = oo for g > g c at T = 0. Hence, the intersection of 
lines for different temperatures gives an estimation of g c . 
By extrapolating the crossing point to T = 0, we obtain 
g c ~ 0.124 with the error 0.001. 



V. SUMMARY 

We have developed an algorithm of CT-QMC for mod- 
els including the spin-boson coupling, i.e., the Bose 
Kondo model and the Bose-Fermi Anderson model. In 
the Bose Kondo model, we have observed the low- 
temperature static susceptibility consisting of the Curie 
term T _1 and the bosonic fluctuating term T~ s . By in- 
cluding the fermionic bath, i.e., in the Bose-Fermi Ander- 
son model, we have identified the quantum phase transi- 



7 



tion between the bosonic fluctuating state and the Kondo 
singlet state. At the critical point, the quasiparticle en- 
ergy scale and the effective moment vanish from each 
side of g. On the other hand, the energy scale of the 
bosonic fluctuation is not affected by the fermionic bath. 
As a result, the static susceptibility exhibits the power- 
law singularity x °< T~ s as the leading term at the critical 
point. 

In this paper, we have restricted ourselves to U = oo in 
the simulation. However, the formalism for the partition 
function, or the weight W(t,/j,), holds also for U < oo, 
so that only the update procedure should be modified 



to take account of the doubly occupied state. One can 
also apply the present framework to the Kondo limit, i.e., 
Bose-Fermi Kondo model. In this case, the algorithm for 
the Kondo model (CT-J) is available. 0, 9 

The present method can apply to lattice models such 
as the Hciscnberg model and the t-J model in terms of 
the extended DMFT equation. An application to lattice 
models as well as a detailed investigation of the impurity 
models are left for future issues. 

The author is supported by JSPS Postdoctoral Fellow- 
ships for Research Abroad. 
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